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Electron-hole generation and recombination rates for plasmon emission and absorption in graphene 
are presented. The recombination times of carriers due to plasmon emission have been found to 
be in the tens of femtoseconds to hundreds of picoseconds range. The recombination times depend 
sensitively on the carrier energy, carrier density, temperature, and the plasmon dispersion. Carriers 
near the Dirac point are found to have much longer lifetimes compared to carriers at higher energies. 
Plasmons in a graphene layer on a polar substrate hybridize with the surface optical phonons and 
this hybridization modifies the plasmon dispersion. We also present generation and recombination 
rates of carriers due to plasmon emission and absorption in graphene layers on polar substrates. 
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The high carrier mobility and the large optical absorp- 
tion in graphene have opened up unique opportunities for 
this material in electronics and optoelectronics 1 - - — . The 
performance of graphene in many of these applications 
depends on the electron-hole generation and recombina- 
tion rates in graphene. It is therefore important to under- 
stand the mechanisms that are responsible for electron- 
hole generation and recombination in graphene and the 
associated time scales. Previously, carrier generation and 
recombination rates in graphene due to Auger scatter- 
ing and impact ionization and due to optical phonon 
emission and absorption have been reporte d 21 ' 22 . The 
strong interaction between electrons/holes and plasmons 
in graphene has been used to explain observed features 
in the angled resolved photoemission (ARPES) data 1 ^— . 
In this paper, we present electron-hole generation and re- 
combination rates due to plasmon emission and absorp- 
tion. Our results show that the recombination times of 
carriers due to plasmon emission are in the tens of fem- 
toseconds to hundreds of picoseconds range and depend 
sensitively on the carrier energy, carrier density, temper- 
ature, and the plasmon dispersion. The available phase 
space for plasmon emission is restricted because of en- 
ergy and momentum conservation requirements and also 
because of Pauli's exclusion principle and carriers near 
the Dirac point have much longer lifetimes compared 
to carriers at higher energies. Plasmons in a graphene 
layer on a polar substrate (Figure ((Hd)) hybridize with 
the surface optical phonons and this hybridization splits 
the plasmon dispersion into two branche o 23 ' 24 . We also 
present generation and recombination rates of carriers 
due to plasmon emission and absorption in graphene lay- 
ers on polar substrates. The results presented here, in the 
light of the previous studie o 21 ' 22 , indicate that plasmon 
emission is the dominant mechanism for carrier recombi- 
nation in graphene. Our results are expected to be useful 
in interpreting experimental observations in ultrafast op- 
tical studies 2 ^— and in understanding the operation of 
graphene based optoelectronic devices 9 - ' 13 ' 38 ' 39 . 



We hrst consider a graphene sheet in the plane z = 
sandwiched by media with free-space permittivity (Fig- 
ure (UK)), as in the case of suspended graphene^ 3 . The 
electron energy dispersion are given by E s (k) = shvk, 
where s equals +1 and -1 for conduction and valence 
bands, respectively. The dispersion for the plasmons in 
given by e(q, ui) = 0, whe re 7 ' 34 ' 35 , 
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The electron-hole propagator H(q, ui) is 7 - ' 34 ' 35 
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The matrix element between the Bloch functions in the 
above expression equal o 7 ' 34 ' 35 , 
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Here, 9 is the angle between k and q. To calculate the re- 
combination and generation rates, we consider a plasmon 
wave with the electric field given by, 
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E(r, z,t) = ~(q ± iz)E a e ' 1 ;l < 
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The transition rate for an electron in the conduction band 
to go into the valence band via stimulated emission of a 
plasmon of wavevector q is given by the Fermi's Golden 
Rule, 
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The energy density W of the plasmon wave has contribu- 
tions from both the field as well as the kinetic energy of 
the carriers. Assuming no plasmon dissipation, the total 
energy density can be found from the complex electro- 
magnetic energy theorem^, 
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Since the conductivity a(q, lu) is related to the dielectric 
constant e(q,ui) as, 



e(q, u>) = l + i 
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the expression for the energy density W becomes, 
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FIG. 1: (a) A suspended graphene sheet, (b) A graphene 
sheet on a polar substrate, (c) Electron-hole recombination 
via plasmon emission in p-doped graphene (p » n). (d) 
Electron-hole recombination via plasmon emission in photoex- 
cited graphene (n = p) . 



W must also equal n(q)huj(q)/ A, where n(q) is the num- 
ber of plasmons in the mode q and A is the area of the 
crystal. Therefore, using (|S|) and ©, the lifetime of an 
electron in the conduction band due to both stimulated 
and spontaneous emission into all plasmon modes be- 
comes, 
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The recombination and generation rates, R and G (units: 
# / cm 2 -s) due to plasmon emission and absorption can be 
written as, 
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In thermal equilibrium, the plasmon number n(q) equals 
the Bose factor (exp(huj(q) / KT) — Similar results 

can be obtained starting from the self-energy of an elec- 
tron in the conduction band. Assuming thermal equilib- 
rium, and using the imaginary-time Green's function ap- 
proach, the relevant contribution to the conduction band 
electron self-energy can be written as^l, 
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where G(k~q 1 uj n — iv m ) is the valence band Green's func- 
tion. The summation over the Matsubara frequencies can 
be performed by first isolating the pole coming from the 
zero of e(q, iv m ) in the denominator at the plasmon fre- 
quency. Finally, if one excludes from (1121) contributions 
coming from plasmon absorption processes and then cal- 
culates the lifetime of the conduction electron using the 
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then the result obtained is identical to the one given ear- 
lier in ©. It should be mentioned here that focusing 
on the collective excitation pole coming from the zero of 
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FIG. 2: The plasmon dispersion (solid) in a p-doped graphene 
is shown for different hole densities. The dashed curve repre- 
sents hvq. T = 30K. 
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FIG. 3: The calculated spontaneous emission lifetime of an 
electron in the conduction band is plotted as a function of the 
electron energy for different hole densities and temperatures 
in p-doped graphene. 



e(q,iu m ) allows one to calculate the interband scattering 
rate due to electron-plasmon interaction. Contributions 
from other processes, such as Auger scattering and im- 
pact ionization^!, are therefore excluded. From the re- 
sults obtained above it follows that the electron-plasmon 
interaction in graphene can be approximately described 
by the following Hamiltonian in the second quantized 
form, 
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where, b$ and c s r are the plasmon and the electron 
destruction operators, respectively, a stands for different 
spins and valleys, and the coupling constant M g g , g _, is 
given by, 
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Here, 9 is the angle between k and q and A is the area of 
the graphene crystal. 



A. Results and Discussion 

The plasmon dispersion is first found numerically us- 
ing the expression for H(q, u) in ([2]). The recombina- 
tion and generation rates and lifetimes are then calcu- 
lated using (|9|), (fl"0|) . and (fTTj) . The dominant contribu- 
tion to the propagator in @ comes from the intraband 
part. The interband part modifies the plasmon disper- 
sion slightly and also imparts an imaginary part to the 



plasmon frequency. If the interband contribution is ig- 
nored the error in the calculated plasmon frequency has 
been found to be generally small (less than 10%) for small 
plasmon wavevectors. This small error comes with the 
enormous simplicity of having to find zeroes of e(q, lo) on 
only the real frequency axis and therefore this approach 
has been adopted in numerical simulations. The results 
we present are also not self-consistent in the sense that 
the quasiparticle density of states have been assumed 
to be that of the non-interacting electron system. It is 
known that electron-plasmon interaction can modify the 
quasiparticle density of states and generate plasmaron 
bands 17 . However, the modification of the quasiparticle 
density of states is expected to be small for the elec- 
tron and hole densities considered in this paper. The 
average recombination and generation times are defined 
as, r^ 1 = R/min(n,p), and r^ 1 = G/mm(n,p), respec- 
tively, where n and p are the electron and hole densities. 

Figure ^ shows the plasmon dispersion in p-doped 
graphene for different hole densities. Figure ([3]) shows 
the calculated lifetime of an electron in the conduction 
band due to spontaneous plasmon emission as a function 
of the electron energy for different hole densities and tem- 
peratures in p-doped graphene (Figure Q}:)). Note that 
the electron-hole symmetry in graphene implies that the 
hole lifetimes in n-doped graphene would be identical. 
Figure ([3]) shows that conduction electrons with energies 
near the Dirac point have much longer lifetimes com- 
pared to the electrons at higher energies. This trend can 
be understood as follows. Energy and momentum con- 
servation require that an electron with wavevector k can 
emit a plasmon with wavevector q only if, 
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The shaded region in Figure (|4|) shows the allowed val- 
ues of the wavevector of the emitted plasmon as a func- 
tion of the electron energy assuming p = 10 12 cm -2 
and T — 30K. As the electron energy becomes smaller 
than ~60 meV, the allowed phase space for the plas- 
mon wavevectors shrinks significantly for small wavevec- 
tors. The probabilities of emission of plasmons of large 
wavevectors are small because the energies of such plas- 
mons are large and the resulting final electron states deep 
inside the valence band are already occupied by valence 
electrons. In addition, numerical simulations show that 

3t^de(q,Lj)/du:\ u> ^\ in (j9|) becomes large whenw(g) ap- 
proaches qv, which happens for very large wavevectors 
and this also reduces the probability of emission of large 
wavevector plasmons. Figure ([3]) shows that the electron 
spontaneous emission lifetimes can range from values as 
small as 10 fs to hundreds of picoseconds. Figure J5]) 
shows the average minority carrier (electron) recombina- 
tion time tr plotted as a function of the minority carrier 
density for different temperatures in p-doped graphene 
(p = 10 12 cm -2 ). As expected from the results in Figure 
([3]), the average recombination time decreases with the 
increase in the temperature because the minority carrier 
distribution spreads to higher energies. 

In many optical studies^— and in graphene based op- 
toelectronic devices ^ 13 ' 38 ! 39 , photoexcitation followed by 
rapid thermalization results in a an equal number of ther- 
mally distributed electrons and holes in an otherwise near 
intrinsic graphene layer (Figure (jTJi) ) . It is therefore im- 
portant to understand the recombination and generation 
times in such situations. Figure (j6|) shows the recom- 
bination times tr (solid) and the generation times tg 
(dashed) plotted as a function of the electron and hole 
density (assumed to be equal) for different temperatures. 
The number of plasmons n(q) in different modes is as- 
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FIG. 4: The shaded region shows the allowed values of the 
spontaneously emitted plasmon wavevectors as a function 
of the electron energy in the conduction band for p-doped 
graphene. p = 10 12 cm" 2 and T = 30K 
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FIG. 5: The average minority carrier (electron) recombination 
time tr is plotted as a function of the minority carrier den- 
sity for different temperatures in p-doped graphene (p — 10 12 
cm -2 ). The arrow indicates curves for increasing values of 
the temperature (T = 30, 100, 300, 900K). 
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FIG. 6: The recombination times tr (solid) and the genera- 
tion times tg (dashed) are plotted as a function of the electron 
and hole density (assumed to be equal) for different tempera- 
tures. The arrows indicate curves for increasing values of the 
temperature (T = 30, 100, 300, 900K). 



sumed to be given by the Bose factor. This assumption 
may not be valid in a non-equilibrium situation imme- 
diately following photoexcitation. Figure ([7]) shows the 
recombination times tr (solid) and the generation times 
tg (dashed) plotted as a function of the temperature for 
different electron and hole densities. Figures (JSJ and ([7J 
show that the recombination times can be much smaller 
than a picosecond for carrier densities larger than 10 11 
cm -2 at all temperatures. Figures ((6j) and (JJJ show that 
the generation times can also be very short and this im- 
plies that carrier generation cannot be ignored in experi- 
ments where a hot carrier distribution is created via pho- 
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FIG. 7: The recombination times tr (solid) and the genera- 
tion times tg (dashed) are plotted as a function of the tem- 
perature for different electron and hole densities (assumed to 
be equal). The arrows indicate curves for increasing values of 
the carrier density (n =p = 10 9 , 10 10 , 10 11 , 10 12 cm" 2 ). 
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FIG. 8: The dispersion of the coupled plasmon-phonon mode 
(solid) in a p-doped graphene sheet on a SiC substrate is 
shown for different hole densities. The dispersion splits into 
two branches. The dashed curve represents Tivq. T — 30K. 
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III. THEORETICAL MODEL FOR GRAPHENE 
ON POLAR SUBSTRATES 

The results presented above suggest that it might 
be possible to alter the plasmon- assisted recombina- 
tion and generation rates in graphene by altering the 
dielectric environment thereby modifying the strength 
of the Coulomb interaction^. Specifically, a substrate 
with a large dielectric constant could potentially reduce 
the recombination and generation rates. However, po- 
lar materials with large dielectric constants have sur- 
face optical phonon modes that couple strongly with the 
graphene plasmons 2 ^. To study this further, we consider 
a graphene sheet at a distance d away from a polar sub- 
strate (Figure (lib))- The dielectric constant of the sub- 
strate is assumed to be given by the expression, 



ib (u) = e sub (oo) 



LO TO J 
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Here, e sub (0)/e sub (oo) = lu 2 lo /lo^ . The surface opti- 
cal phonon frequency to so is obtained by setting e sub (uj) 
equal to -1, and equals, 



' e SMfc (0) + 1 
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The dielectric constant e(q, co) of the graphene sheet can 
be found by placing a test charge in the sheet and finding 
the resulting potential. The result is, 
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The dispersion of the coupled plasmon-phonon longitudi- 
nal mode can be found as before by setting e(q, ui) equal 
to zero. Now one finds two longitudinal collective modes. 
In the q — > limit, the lower frequency mode is plasmon- 
like with ui(q) — > as q — > 0, and the higher frequency 
mode is phonon-like with u)(q) — > ojso as q — > 0. For large 
wavevectors, the lower frequency mode disappears into 
the electron-hole continuum while the higher frequency 
mode becomes plasmon- like with tu(q) — > qv as q — > oo. 
As an example, we consider the technologically relevant 
case of a graphene layer on a Silicon Carbide (SiC) sub- 
strate ^ 23 ' 24 . The values of different parameters are as 
follows: d — 5 Angstroms, Uujlo = 120 meV, Tiujto = 98 
meV, and = 6.5^2. These give Kujso ~ H7 meV. 
Figure ((5J) shows the dispersions of the coupled plasmon- 
phonon modes for a p-doped graphene sheet on a SiC 
substrate for different hole densities. Comparing Figures 
© and flSJ), it can be seen that plasmon-phonon cou- 
pling significantly modifies the dispersion and this has 
recently been verified experimentally 23 . The recombi- 
nation and generation rates can be obtained using the 
same expressions as those given in ([§]), (JTUJ), and pT|) 
with the exception that contributions from both branches 
of the dispersion must be included. Therefore, surface 
optical phonons of the polar substrate provide an addi- 
tional channel for carrier recombination and generation. 
It should be mentioned here that large wavevector surface 
optical phonon modes can also cause intervalley recombi- 
nation and generation processes 2 ^. However, the square 
of the coupling matrix element between the surface opti- 
cal phonons and the carriers is proportional to^ 2 -, 
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FIG. 9: The calculated spontaneous emission lifetime of an 
electron in the conduction band is plotted as a function of the 
electron energy for different hole densities and temperatures 
in p-doped graphene on a SiC substrate. 




FIG. 11: The recombination times tr (solid) and the gen- 
eration times tg (dashed) are plotted as a function of the 
electron and hole density (assumed to be equal) for different 
temperatures for a graphene sheet on a SiC substrate. The ar- 
rows indicate curves for increasing values of the temperature 
(T = 30,100,300,900K). 
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FIG. 10: The average minority carrier (electron) recombina- 
tion time tr is plotted as a function of the minority carrier 
density for different temperatures in p-doped graphene (p = 
10 12 cm -2 ) on a SiC substrate. The arrow indicates curves for 
increasing values of the temperature (T = 30, 100, 300, 900K). 



and becomes small for the large wavevectors needed for 
the intervalley transitions in graphene. Therefore, inter- 
valley processes will be ignored here. 



A. Results and Discussion 

Figure © shows the calculated lifetime of an electron 
in the conduction band due to spontaneous emission as a 
function of the electron energy for different hole densities 
and temperatures for a p-doped graphene sheet on a SiC 
substrate. Figure © displays the same general trends as 



does Figure @ in the case of a suspended graphene sheet. 
However, lifetimes are shorter for the low energy electrons 
in the case of graphene on SiC. For electrons with ener- 
gies near the Dirac point, recombination is entirely due 
to the lower frequency branch of the dispersion which 
facilitates interband transitions better than the plasmon 
dispersion in suspended graphene. The sharp peaks seen 
in Figure © occur when the lifetimes due to the lower 
frequency branch of the dispersion are becoming longer 
with the electron energy while lifetimes due to the up- 
per frequency branch are becoming shorter. As in the 
suspended graphene case, the spontaneous emission life- 
times can range from tens of femtoseconds to hundreds 
of picoseconds. Figure (|10p shows the minority carrier 
(electron) recombination time tr plotted as a function of 
the minority carrier density for different temperatures in 
a p-doped graphene (p = 10 12 cm~ 2 ) on a SiC substrate. 
Compared to the suspended graphene case (Figured} , the 
recombination times for graphene on SiC are shorter for 
small minority carrier densities. Next, we consider the 
case when the electron and hole densities are the same 
(as is the situation in photoexcitation experiments). Fig- 
ure (|lip shows the recombination times tr (solid) and 
the generation times tq (dashed) plotted as a function 
of the electron and hole density (assumed to be equal) 
for different temperatures for a graphene sheet on a SiC 
substrate. The recombination and generation times in 
graphene on SiC are generally of the same order as in 
the case of suspended graphene discussed earlier. The 
role of the higher dielectric constant of the SiC substrate 
in reducing plasmon-assisted recombination and genera- 
tion rates, compared to suspended graphene, is compen- 
sated by the presence of surface optical phonons which 
not only modify the plasmon dispersion but also provide 
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an additional channel for recombination and generation. 

IV. CONCLUSION 

In this paper we have presented electron-hole re- 
combination and generation times due to spontaneous 
and stimulated emission and absorption of plasmons in 
graphene. Our results indicate that plasmon assisted re- 
combination times in graphene can vary over a wide range 
of values ranging from tens of femtoseconds to hundreds 
of picoseconds. In many proposed and demonstrated op- 
toelectronic devices ^ 1 ^ 13-1 ^ 38 ' 39 , the plasmon-assisted 
recombination and generation rates could be fast enough 
to significantly impact device performance. 
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